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We derive a simple formula for the finite-temperature shifts 
of the collective excitations of a Bose-Einstein condensate. To 
test the validity of our treatment we apply it to explain the 
anomalous behavior of the "m = 2" excitation frequency in 
the recent JILA experiments (D. S. Jin et al, PRL, 78 (1997) 
764), where this frequency does not approach its "noninter- 
acting" value at u = 2u> p but, instead, decreases with tem- 
perature. It turns out that the effect is mainly governed by 
a resonance between the condensate excitation frequency and 
one of the quantum eigenfrequencies of the thermal cloud. 
Comparison of our predictions with the JILA experimental 
data shows a good agreement. 

The recent discovery of the Bose-Einstein condensa- 
tion in dilute atomic gases jjj has incurred a tremen- 
dous burst of experimental and theoretical research in 
this area. Whereas zero-temperature properties of Bose 
condensates seem to be well understood by now, theory of 
the finite-temperature behavior of the system is far from 
being complete. One of the most intriguing effects unex- 
plained so far is the shifts in the spectrum of collective 
excitations at finite temperatures. A possibility of direct 
comparison of the theoretical predictions in this area with 
the experimental data provides us with a unique op- 
portunity to test the numerous theoretical assumptions 
used to describe Bose gases in general. Among such open 
issues are validity of the mean-field description, relevance 
of the anomalous density of the non-condensed particles, 
role of the many-body and trap-induced modifications of 
the interaction strength etc. 

The first attempt to approach the problem of the finite- 
temperature shifts was made in the works [^,^) where the 
collective excitation frequencies were associated with the 
elementary excitations spectrum of a stationary Popov 
Hamiltonian, where the thermal cloud was considered as 
being static. In a subsequent paper S the model was im- 
proved by taking into account the many-body reduction 
of the interaction strength fflj. The next important step 
forward involved the dynamical response of the thermal 
cloud |^||. Furthermore in the Ref. |10| the frequency 
shifts were calculated using Gaussian ansatz for both 
the condensate and the thermal cloud. Hydrodynamic 
treatment of the problem is presented in Whereas 
all the works listed above were based on the mean-field 
treatment authors of the paper jl2| develop a perturba- 
tion theory to derive expressions for both frequencies and 
damping rates. 

In our work we first derive a general formula for the 



finite temperature frequency shifts of the collective ex- 
citations of a Bose condensate: the expression we ob- 
tain is independent of model chosen to describe the ther- 
mal cloud. Furthermore we notice that the recent finite- 
temperature excitation experiments in JILA [Q give us 
an excellent opportunity to test our findings: using the 
fact that in the case of the "m=2" excitation a strong res- 
onance between the collective excitation of interest and 
one of the quantum eigenfrequencies of the thermal cloud 
takes place, we restrict our description of the dynamical 
response of the thermal cloud to only two modes and 
obtain a simple closed set of algebraic equations, from 
which we deduce the values of the frequency shifts at dif- 
ferent temperatures. We compare our results with the 
JILA experimental data. 

Let us consider a time-dependent Gross-Pitaevskii-like 
equation for the mean cf> = (i/)) of the quantum field if) of 
a collection of trapped bosons below the Bose-Einstein 
transition temperature: 



+ g(n c + 2n) }<f> + grh<f>* ; (1) 





here and below n c = N c \(j>\ 2 and m c = N c (j> 2 are the 
condensate normal and "anomalous" densities, iV c is the 
number of particles in the condensate, n = (Stp^St/j) and 
fh = {S'l/jS'tp) are the normal and "anomalous" densities 
of the excited particles expressed through the variance of 
the field 5tp = ip — (if)) , the condensate wave function is 
supposed to be normalized to unity at the initial instant 
of time, coupling constant g is defined as g — 4irh 2 a/M, 
a is the s-wave atomic scattering length, and M is the 
atomic mass. Note that if the thermal densities h and rh 
are known the above equation is exact |i~5| ] . 

We will be interested in oscillating solutions of (|l|) 
whose frequency is close to the frequency Wk moda of the 
one of the zero-temperature condensate excitations. We 
will use the following ansatz 

m (2) 

+ £)(/) (_) e _l(a ' k modc+ £ ) t + D*0(+) e + i ("k modc +e)t| 

n(t) = n (0) + n(-) e - l (^„,odc+ £ ) 4 + n(+) e +^* modc + £ )t (3) 
m(t) (4) 
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to describe the motion of the condensate and thermal 
cloud. Here e is the frequency shift of interest, // is the 
chemical potential, and D is the amplitude of the conden- 
sate oscillations. We will suppose also that the thermal 
densities n and rh are known functions of the amplitude 
D: 



n(D) 



m(D) 



(5) 



Furthermore we will assume that the oscillations of the 
thermal density are "induced" by the condensate oscilla- 
tions and, therefore, the (+)- and (— )-components of the 
thermal density go to zero as the condensate oscillation 
amplitude D goes to zero: h^~' + \D = 0) = jfi^' + \D = 
0) = 0. 

There are two independent small parameters present in 
the problem: the ratio between the non-condensed and 
condensed densities n/n c and the amplitude D. Let us 
now expand </>, /x, and e into a double series with respect 
to both h/n c and D, expand h and m into a simple series 
with respect to D, and then insert this expansions in 
the condensate equation ([!]). To zeroth order in h/n c 
and for the infinitely small amplitudes of the condensate 
oscillations, the stationary part of the condensate wave 
function is given by the ground state solution <!> of the 
Gross-Pitaevskii equation 



2M + VtTap 



■gN c \<s>\ 



<f> = Xi <I> 



(0) 



where the chemical potential /xo is chosen in such a way 
that the condensate wave function is properly normal- 
ized: ($|$) = 1. The oscillating part of the condensate 
wave function 



V (</>(+))* 



(7) 



appears only in the first order in D (and zeroth order in 
h/n c ): it is given by 
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yt 



(8) 



where |W / k modc y= (£^k modo , Vk modo ) is the mode func- 
tion of the collective excitation we work with. This 
mode function is represented by an eigenstate of the 
Bogoliubov-de-Gennes Liouvillian 



C = 



Mn 



-MX -Ln 



(9) 



L ° = 2M + Vtiap + 2gNc ^ 2 ~ fl ° 
M = gN c t> 2 

of an energy -Ek mode = fiwi tmodi , . The eigenstate 
|Wk modo y is supposed to be "normalized" to unity: 
~^Wk mode \Wk mod< ,y= 1 , where the bilinear form -<-|«>- 
is defined as 



^w p \w q y= J d 3 r {u;u q - v;v q ) 



(10) 



Let us now consider the terms of the first order in 
h/n c and zeroth order in D in the lhs and rhs of the 
condensate equation (Q) . To this order of approximation 
we can calculate the thermal-cloud-induced correction to 
the chemical potential A/x = \i — /iq, and it reads 



A/x =^ a \A£ (< $ =0 \a 3 $y 



(11) 



Note that this correction is given by an off-diagonal ma- 
trix element of the stationary part of the thermal correc- 
tion to the unperturbed Liouvillian 



a3 q] 



-g(m<®)* -2g(n<®) 



(12) 



calculated between the zero-energy "condensate" eigen- 
mode |$ >-= — $*) and a state |$„ y= 



UN, 



(VN2cr 3 \^ y) introduced in @||| (so-called 



"lost eigenvector": it is the only eigenvector of the op- 
erator (/Co) 2 not present in the set of the eigenstates of 

Cq itself). The matrix 173=' ' 



( s 1 , is the third Pauli 
matrix. 

To the same order in n/n c and D one can show that the 
correction to the condensate wave function A<1> = (f>^ — $ 
is given by a solution of an inhomogeneous equation 



£ |A^= -(A£ ( 2 =0 - A/x) \a 3 <S>y 
and is a subject of an additional constraint 
^$|A<£>^= . 



(13) 



(14) 



Here the vector |A$^ is defined as | A$^= (A<I>, A$*). 

Consider now the terms in the equation (|l|) which are 
linear with both n/n c and D. To this order the frequency 
shift we are interested in is given by 

Tie = ^Wk mod JA4°) + A£ ( 2 =0 - AfX(x 3 \Wk moda y (15) 
+ ^ kmodc |A A £^+)|a 3 ^ 

where the condensate-induced correction to the Liouvil- 
lian (|)|) and the oscillating part of the thermal correction 
to it read 

AC C = (16) 

2g[$*A$ + A$*$] 2g$A$ 

-2g$*A$* -2 ff [$*A$ + A$*$] 

-2j(8W) 4 



A £(-/+) 



(17) 



respectively. The formulas ( Ofl^ , ff^ ) defining the finite- 
temperature frequency shift of the condensate oscillation 
frequency is the central result of our paper. 
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In the derivation above we have used a conjecture 
of completeness of the basis formed by the vectors 
{ |$ a ^; |W k ^; } (see @), where the 

mode functions |Wk>*~ are the positive energy eigenstates 
of the unperturbed Liouvillian: 



C \W k y= E k \W k y 



(18) 
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and o\ = ^ ^ q J is the first Pauli matrix. Note that to 

the first order in h/n c the thermal-cloud-induced correc- 
tions to the mode energies read 



;(o) 

D=0 



An<j 3 \W k >- 



(19) 



To define the frequency shift ([L5|) in an excitation ex- 
periment one should know the stationary part C 2_ of 
the non-condensed corrections to the condensate Liouvil- 



lian and the "dynamical polarizability" AC' D =o ■ 



In 



general such a calculation would require a detailed knowl- 
edge of the dynamical response of all the thermal modes 
and it is a complicated numerical problem. However in 
some cases the calculation can be simplified: that is, for 
example, the case of an accidental degeneracy between 
the condensate oscillation frequency Wk modo and one of 
the eigenfrequencies of the thermal cloud Wk 2 — ^ki '■ 



^k n 



w k 2 



Wki 



(20) 



In this case the time evolution of the cloud gets domi- 
nated by the two resonant modes (1 and 2) and the ther- 
mal cloud can be treated in a kind of resonant approxima- 
tion, where only two modes evolve with time. The latter 
assumption is justified for the case of the "to = 2" ex- 
citation experiment at JILA where the frequency of the 
mode of interest Wk modc = uJ n =a,m=+2 ~ 1-45 u> p turns 
out to be close to the differential frequency — w kl = 

4^2 m=+i — w n=0,m=-i ~ 1.65 (mode classification 
corresponds to one suggested in |16[ ). 

At this point we have to chose a specific model to de- 
scribe the evolution of the thermal cloud: in what fol- 
lows we will treat the thermal particles using a time- 
dependent Hartree-Fock-Bogoliubov approximation Q. 
Under this approximation the evolution of the variance 
Sifj — xp — (ip) of the bosonic field will be governed by a 
mean-field equation 



4^ 



2M 
g{m c 



^trap + 2g(n c + n) } Sip 
- to) 5$' . 



(21) 



Furthermore we will suppose that prior to the excitation 
the system was in a thermal equilibrium. 

Following the "resonant approximation" program de- 
scribed above one can show that in the resonant case 
the evolution of the thermal cloud is governed by a sim- 
ple equation for an effective two-level atom in a resonant 
field of a Rabi frequency 



hn Rahi /2 = 2gN c ^M/ k2 |{(a 3 $) ® W kmodc }| W kl > D 

+0(D 2 ) (22) 



and detuning 
US = %e- [{E k2 + A£k 2 



(£ kl +A£ kl )-£ k „J 

+0(D) , (23) 



where the energy shifts AE k are given by (19) and the 
tensor product {• ® •} is defined as 



{W p ® W q } 



(24) 



[u;u q + v;v q ] 



-u;v q 



v;u q 



-p;u q 



v*v q ] 



We can identify two distinct classes of small collective 
excitations of the system: for the excitations of the first 
class (|) the stationary part of the thermal cloud is the 
same as in thermodynamic equilibrium prior to the be- 
ginning of the excitation process; for the second class (t) 
the populations of the two resonant modes are inverted 
with respect to the thermal equilibrium. 

The thermal Liouvillian components (|l^, [Tt]) we are 
interested in now read 

AC %— — A£ cqu iiib. (25) 
+2g(N kl - N U2 ) ({W k2 ® W k2 } - {W kl ® W kl }) 

'w+V 



— A£ ( -/+ 3 = 2g(N kl - N k2 ){W kl <g> W k J 



d 



x ( W P2i ) , (26) 



where 



, for the (|)-mode 
, for the (t)-mode 



is the population inversion; 



P21 



"^RabiW / 25 



(27) 



(28) 



is the coherence induced between mode 2 and mode 1; 
7V k = {exp(£k/T) — 1} _1 are the Bose-Einstein occupa- 
tion numbers in the thermal equilibrium; A£ oqu iiib. = 
((2#n cqui i ib ., 3TO*cquiiib.), (-gm* quilih , -2gn eqai i ih .)) is 
the thermal cloud induced correction to the Bogoliubov- 
de-Gennes Liouvillian (||) calculated in the state of ther- 
mal equilibrium. Notice that we have preserved the first 
order in n/n c terms in the denominator of the expression 
for the detuning fl25|), though this seems to contradict 
the expansion procedure we have chosen. In the resonant 
case, however, the presence of this terms is well justified 
since in this case the mean-field shifts A£k and the fre- 
quency shift e can approach the "zeroth order" splitting 
E^ 2 - E kl - i? kmodc . 
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FIG. 1. Shift of the frequency of the "m=2" collective ex- 
citation with respect to its zero-temperature value as a func- 
tion of temperature. Lines show the collective modes we have 
found. Down(up) arrow corresponds to the normal(inverted) 
mode. Errorbars show JILA experimental data [2] . As a unit 
of frequency the radial trap frequency ui p has been chosen. 
The inset illustrates the relationship between the energies of 
the modes involved in the model. 



The formula for the frequency shift ( |15| ) along with 
those for the corrections to the Bogoliubov-de-Gennes Li- 
ouvillian ([l6p5|) and dynamical polarizability of the ther- 



mal cloud (p6[) provides us with a closed implicit system 
of algebraic equations from which one can extract the val- 
ues of the finite-temperature resonant corrections to the 
oscillations frequencies of the condensate. At the Fig.[l] 
we show a comparison of our predictions with the exper- 
imental measurement of the "m=2" condensate frequen- 
cies at JILA. We would suppose that in the temperature 
range below T ~ .65 T c the "normal" mode was excited 
in the experiment, whereas above this temperature the 
experimental points correspond to the "inverted" mode. 
The point at T ~ .65 T c can be associated with any of the 
two modes or, probably, with a superposition of them. 

In our calculations the condensate wave- function $, 
"lost vector" $ a , chemical potential \i, and thermal- 
cloud-induced correction to the condensate wave- function 
A$ were described within the Thomas-Fermi approxi- 
mation. The mode functions |Wk mode >-= | W n =o,m=+2'^ , 
|W kl >-= \W n=0 , m =-i>, and \W k2 y= \ W„ =2 , m =+i,[-] >- 
were approximated by the Thomas-Fermi expressions 
presented in | Ea |. For the mode energies we were using 
the finite- N numerical values calculated in |l7]] , as the ef- 
fect we are describing is very sensitive to the actual value 
of the energy splittings. Throughout our calculations we 
where neglecting the thermal equilibrium part A£ equ iiib. 
of the Liouvillian as this type of perturbation has been 
shown to have no significant effect in the frequency shifts 



As we pointed out above, one of our collective modes 
corresponds to the inversion of populations of the ther- 
mal modes 1 and 2. Whether such an inversion can be 
reached in the experiment depends on details of the exci- 
tation procedure. We have checked, however, that the 
3% amplitude variation of the trap frequency used in 
the JILA experiment is a strong saturating perturba- 
tion of the transition between the {n = 0,m = — 1} 
and {n = 2,m = +1, [— ]} thermal modes, and, there- 
fore, the "inverted" mode could, indeed, be excited dur- 
ing the excitation stage. The answer of the question on 
which particular mode will be excited at a given tem- 
perature would involve a time-dependent analysis of the 
whole experimental sequence: it is a subject of our future 
research. 

Two conclusions can be drawn from our work: (a) At 
finite temperatures the spectrum of the collective exci- 
tations of the condensate is essentially different from the 
spectrum of the elementary excitations, in contrary to the 
zero-temperature case. Comparison of the expressions 
( jl5| ) and ( p^ ) for the finite-temperature corrections to 
the former and latter shows that the collective excitation 
shift contains an extra term proportional to the "dynami- 
cal polarizability" of the thermal cloud; (b) Quantization 
of motion of the thermal quasi-particlcs is important and 
it can lead to significant effects, at least in the resonant 
case described in our paper. 

Note also, that in the off-resonant limit our results are 
consistent with the conclusions of the work |8| . 
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